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Abstract 

We study the evolution of spherically symmetric radiating fluid dis- 
tributions using the effective variables method, implemented ab initio in 
Schwarzschild coordinates. To illustrate the procedure and to establish 
some comparison with the original method, we integrate numerically the 
set of equations at the surface for two different models. The first model 
is derived from the Schwarzschild interior solution. The second model is 
inspired in the Tolman VI solution. 

1 INTRODUCTION 

In 1980, Herrera, Jimenez and Ruggeri proposed a scminumerical method, here- 
after referred to as the effective variables method (EVM), that can be used to 
obtain nonstatic models from static solutions. This method divides the space- 
time in two spatial regions. The outer region is described by the Vaidya solu- 
tion and the spacetime metric in the interior is obtained by solving the Einstein 
field equations. Further, proper boundary conditions are imposed in order to 
guarantee a smooth matching of both solutions. The EVM has been used ex- 
tensively to study astrophysical scenarios in radiative coordinates |l], ^]. The 
method assumes that the called effective variables p and p, which depend also 
on the time-like coordinate, have the same radial dependence as the corre- 
sponding static physical variables (energy density and pressure) obtained from 
a static interior solution of the Einstein equations. The rationale behind such 
an assumption is the fact that the effective variables reduce to their physical 
counterparts in the static limit. This approach can be justified by means of 
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the characteristic times for different processes which take place in the collapse 
scenario ||, 0, |1 . If the hydrostatic time scale, Th y dr. ~ 1/ VGp, is much shorter 
than the Kelvin-Helmholtz time scale tkh, then in a first approximation the 
inertial terms in the equation of motion can be ignored. Consequently, it seems 
reasonable to assume, in this approximation, that the radial dependence of the 
physical variables is the same as in the static solution. However, a better approx- 
imation is obtained by assuming that the effective variables, not the physical 
ones, have the same radial dependence as the corresponding physical variables 
of the static situation . 

If the EVM is truly general, an implementation in Schwarzschild coordinates 
could be more interesting for astrophysicists since these are the type of coordi- 
nates commonly used by them; this is successfully accomplished in this paper. 
The idea that we can always construct dynamical solutions from static ones 
seems a general method. Besides, as we will show, the EVM in Schwarzschild 
coordinates introduces higher dynamic corrections, by means of the velocity re- 
spect to a minkowskian observer, than those obtained in the Bondi coordinates 
treatment. 

In this paper we model two simple but physically meaningful scenarios. The 
first is obtained from a static solution for an incompressible fluid. Once the 
sphere departs from static equilibrium by the emission of energy, it slowly re- 
covers the initial state (staticity), unless some energy can be reabsorbed at the 
surface to reach a constant radius in short time. In fact we show that even in 
Bondi coordinates the same result holds if staticity at the surface is enforced. 
The second model corresponds to a highly compressed gas of fermions, Tolman's 
solution VI H , and leads us to an exploding sphere as has been reported |], [?J . 
Additionally, these seminumerical models could serve as useful test beds for the 
numerical relativity methods being developed to match Cauchy and character- 
istic codes 

The paper is organized as follows. In section 2 we write the field equa- 
tions, for the inner region, in Schwarzschild coordinates. Describing the exterior 
space-time by means of the Vaidya metric, in section 3, we treat the matching 
conditions and write the equations at the boundary of the distribution of mat- 
ter. In the section 4 we show our two models; finally we conclude in section 
5. 



2 FIELD EQUATIONS 

To write the Einstein field equations, inside the distribution of matter, we 
use the line element in Schwarzschild coordinates 

ds 2 = e u dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 d(/) 2 ) , (1) 

where v — v(t, r) and A = X(t, r), with (i, r, 9, </>) = (0, 1, 2, 3). 

Physical input is obtained by introducing Minkowski coordinates (r, x,y, z) 
by@ 

dr = e u/2 dt, dx = e x/2 dr, dy = rd9, dz = r sin 6d<j>. (2) 
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In these expressions v and A are constants, because they only have local values. 

Next we assume that, for an observer moving relative to these coordinates 
with velocity u> in the radial (x) direction, the space is filled with a fluid of 
density p, pressure p, and unpolarized radiation of energy density e. For this 
comoving observer, the covariant energy tensor in Minkowski coordinates is 



/ p + e 
—e 





(3) 



-e 0^ 
p + e 
p 
p J 

Note that from (0) the velocity of matter in the Schwarzschild coordinates 



is 



dr 
~dt 



= Loe^' 2 . 



(4) 



Now, by means of a Lorentz boost and defining e = e(l +w)/(l — uS) we write 
the field equations in relativistic units (G = c = 1) as follows: 
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where the comma subscript represents partial differentiation with respect to the 
indicated coordinate. 

We have four field equations for four physical variables (p, p, e and a;) and two 
geometrical variables (y and A). Obviously, additional information is required 
to handle the problem consistently. First, however, we discuss the matching 
with the exterior solution and the surface equations that govern the dynamics. 



3 MATCHING CONDITIONS AND SURFACE 
EQUATIONS 



We describe the exterior space-time by the Vaidya metric 
2M(u)' 



ds z 



1 



R 



du 2 + 2du dR - R 2 (dO 2 + sin 2 9 1 



0) 
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where u is a time-like coordinate so that u = constant represents, asymptoti- 
cally, null cones open to the future and R is a null coordinate {gun. = 0). 

The exterior and interior solutions are separated by the surface r = a(t). To 
match both regions on this surface we require the Darmois matching conditions. 
Thus, demanding the continuity of the first fundamental form, we obtain 

^• = 1-^ do) 

and 

Va = ~Xa- (11) 

From now on, the subscript a indicates that the quantity is evaluated at the 
surface. Matching conditions are usually obtained from the continuity of the 
first and second fundamental forms. Here, however, we will use the continuity 
of the independent components of the energy-momentum flow instead of the 
second fundamental form, which have been shown to be equivalent \ \n\ but it 
is simpler to apply in the present case. This last condition guarantees absence 
of singular behaviors on the surface. It is easy to check that 

Pa = 0, (12) 

which expresses the continuity of the radial pressure. 

To write the surface equations we introduce the mass function m by means 

of 

e - A ( r <*> = l-2m(r,t)/r. (13) 
Substituting ( |l3| ) into (||) and (||) we obtain, after some rearrangements, 



dm 
~dt 



- P + e(l -*,)(!--) 



(14) 



This equation shows the energetics across the moving boundary of the fluid 
sphere. Evaluating (|l^) at the surface and using the boundary condition (|l2|), 
the energy loss is given by 

rh a = -47ra 2 e a (l - 2m a /a)(l - u a ). (15) 

Hereafter, a dot over any variable indicates d/dt. The evolution of the boundary 
is governed by equation (||) evaluated at the surface 

a = (1 - 2m a /a)uj a . (16) 

Scaling the total mass m ai the radius a and the time-like coordinate by the 
initial mass m a {t = 0) = m a (0), 

A = a/m„(0), At = m a /m a {0), t/m a [0) -» t, 

it is convenient to define 

F = 1 - 2M/A, (17) 
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= w„. (18) 
Also we define the luminosity as seen by a comoving observer as || 

£ = (47rr 2 e) r=a , (19) 
and the luminosity perceived by an observer at rest at infinity as 

L = —M = FE(1 + ft). (20) 
The function F is related to the boundary redshift z a by 

1 + Za = = F -V2. (21) 



Thus the luminosity as measured by a noncomoving observer located on the 
surface is 

E = L{l + z a f = -— = E(1 + Q), (22) 
r 

where the term (1 + f2) accounts for the boundary Doppler shift. With these 
definitions the surface equations can be written as 

A = FQ, (23) 

p={ l-F)A + 2L_ 
A 

Equations ( p3| ) and ( p4| ) are general within spherical symmetry. We need a 
third surface equation to specify the dynamics completely for any set of initial 
conditions and a given luminosity profile L(t). For this purpose we can use the 
conservation equation T^t — evaluated at the surface. After straightforward 
manipulations the condition Tj: = results in 

(p + p)(47rr 3 p + to) 



P,r + 



r(r — 2m) 



m,u + —^--^r 1 )+-(p-P), (25) 



47rr(r — 2m) \ r — 2m 2 I r 

where the effective variables are defined by 

and 

__ p + puj' 2 

P = T — T + e. (27) 

These effective variables are essentially the same as have been defined by Her- 
rera and colaborators, but now the velocity to introduces a higher dynamics 
correction (quadratic). This fact could be of interest to investigate its effect on 
dissipative processes (like heat flow and viscosity). 

Equation (^5|) is the generalization of the Tolman-Oppenheimer-Volkov equa- 
tion for hydrostatic support in nonstatic radiative situations. Our equation 
leads at the surface to a differential equation for Q if we specify in some way 
the geometrical variables. 
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4 MODELING 



From (g), (|) and @, easily we obtain 

m = [ 4Trr 2 pdr, (28) 
Jo 

r 2(47rr 3 p + m) , 
J a r(r - 2m) 

Thus, m and v are expressed in terms of p and p in the nonstatic case in the 
same way they are in terms of p and p in the static case. These considerations 
suggest the application of the EVM which until now has been exclusively used 
in Bondi coordinates |l], 0]; that is, we assume that the r dependence on p and 
p is the same as on the p static and p static- 

To illustrate the procedure, in what follows we model two simple scenarios 
which correspond to an incompressible fluid and to a highly compressed gas of 
fermions. 



4.1 Schwarzschild— like model 

Consider the well known Schwarzschild interior solution, where the density sat- 
isfies p = constant. Thus, in the EVM we take the effective density as 

P = /(*), (30) 

where / is an arbitrary function of t. Now, with p = constant we can integrate 
equation (25) in the static case and obtain the expression for p which leads us 
to 

*±2l = (1 - ^pr^m, (31) 
p + p 3 

where k is a function of t to be defined from the boundary condition ( |l2| ) , which 
now reads, in terms of the effective variables, as 

p a = p a n 2 + e a (l-n 2 ). (32) 

Thus, @ and @ give 

3(1 — F) 

P = o_„2 ' ( 33 ) 



8na 
s/F_ 

3 1 - xsVf J 



- = gf Xg VF-3^\ (34) 



with 

£ = [l-(l-F)(r/a) 2 ] 1/2 

and 

Xs = 3(tt 2 + 1)(1 - F) + 2E(1 + 0), 
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Using 



tps = {3SI 2 + 1)(1 — F) + 2E(1 + SI). 
and ( p^ ) it is easy to obtain expressions for m and v. 

m = m a {r/af, 

2 

I XsVF - il> 3 £ | 

e 



(35) 
(36) 



{ 2(1 -F) J 

In order to write down explicitely the surface equations for this example, it is 
interesting to note that the left side of ( p5| ) is zero for any value of f(t). Next, 
evaluating (E5h at the surface, we obtain 



+ (^)](™ + f< 1 -^ 
6FEA 6(1 - F)i 2 



.4 



FE-FE 



-l(l-F)FSl-^ s -{l-F)] 



(37) 



This last equation, together with (£3j) and (|24|), constitute the differential system 
for the surface in this example. It is necessary to specify one function of t 
and the initial data. To compare with Ref. jjj we choose L to be a gaussian 
and radiating away 1/10 of the initial mass. Therefore, the system can be 
numerically integrated for the following initial conditions (among others): 

A(0) = 5.0, F(0) = 0.6, SI = 0.0. 

The integration was done up to some t with good behavior in the physical 
variables. 

Feeding back the numerical values of A, F and SI (and their derivatives) in 
( |35"| ) and ( |3^ ) we obtain m and v (and their partial derivatives) for any value of r. 
Thus, functions p, p, dr/dt and e can be monitored for any piece of the material, 
via field equations. We calculated them for the values r/a — 0.2, 0.4, 0.6, 0.8 
and 1.0. It is interesting to observe that once the sphere departs from static 
equilibrium by means of an emission of energy, it will not recover a constant 
radius at all, at least within the integrated interval of time. This behavior is not 
evident in Bondi coordinates, as we show below. Figure 1 shows the evolution of 
the radius A in a logaritmic scale. Figure 2 displays the profiles of the physical 
variables versus the time-like coordinate for the different comoving regions. 

In order to explore the slow recovery of staticity at the surface we force it 
to return to rest quickly after the emission of energy; we prescribe its evolution 
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instead of giving the luminosity L. For such a prescription the surface equations 
change; now we obtain the luminosity profile from a differential equation. We 
choose a radius evolving as 

{Ai -A f )(e-^ + l) 
A(t > - e (t-U)/a+l + A f ' 

where A^ is the initial radius, Af the final radius, td the decay time and a is 
the decay rate. Figure 3 shows the luminosity profile and the prescribed radius. 
Observe that it is necessary to absorb some quantity of energy (less than the 
total emitted) to reach a final constant radius. We confirmed the same behavior 
in Bondi coordinates. It is interesting to note that in Ref. jD the interior 
profiles of the flux radiation (for the Schwarzschild type model) have the same 
qualitative behavior as we show in Figure 3: It is clear the absorption of energy. 



4.2 Tolman VI— like model 

In this subsection we discuss the model obtained from Tolman's solution VI Q, 
[||. Let us take 

~P=% (38) 
. g jl__9Dr\ 

p = -*\i—d7[ (39) 

where g and D are functions of t, which can be determined using ([32]). Thus, 

* 3(1 -.F) 



247rr 2 ' 



(40) 



-_ (l-F) f ^T-9XT(r/a) \ . 
P 247rr 2 \ i> T - Xr(r/a) J ' [ ' 



where 

Xt = (3^ 2 - 1)(1 — F) + 6E(1 + Q), 



ip T = 3(f> 2 - 3)(1 - F) + 6E(l + Q). 
Using (|8|) and ( p9| ) we obtain 

m = m a r/a, (42) 

{/ \ 2 n 4(1-F)/3F 

In this model, the LHS of (^B|) results in 

Ft(t\ = (! ~ F ) f (jh - 9xt)xt _ , 9xr \ 
n{) l2na^ T -XT)\ 2(^t-Xt) VT 2 J + 

' [(1- F){l + n 2 ) + 2E(l + n)} 2 , (44) 



167ra 3 F 
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which lets us write the third equation at the surface as 



n = 



2 f F 2 A 




|f--[(1- j F)(0 2 + 1) 



+2E(1 + fi)]} - {£ 




(45) 



Again the system can be numerically integrated for a reasonable set of initial 
conditions (with L being a gaussian as in the Schwarzschild-like model, radiating 
away 1/100 of the initial mass), as for example, 



Figure 4 shows the evolution of the radius of the sphere. Initially the fluid 
sphere collapses, but later it bounces. Figure 5 gives the evolution of the matter 
variables for different regions. Note, for the matter velocity profiles, that some 
inner zones continue contracting after the bounce of the outermost ones. 



We have sought the dependence of the EVM upon the Bondi coordinates. The 
idea that we can always construct dynamical solutions from static ones seems a 
general method. In this paper we accomplish successfully such a construction, 
at least, for the Schwarzschild coordinates. Some research is in progress for more 
realistic luminosity profiles and dissipative transport mechanisms, considering 
extended thermodynamics under the time relaxation approximation. We con- 
sidered in this paper two simple and idealized models, not deprived of physical 
meaning at all. We may hope that they contain some of the essential features 
of gravitational collapse inasmuch as we have fed the models with some obser- 
vational data (initial velocity and total mass radiated.) Also, we can hope that 
our toy models could serve as test beds for the numerical relativity methods 
and codes. 
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A(0) = 6.67, F(0) = 0.70, ft = -0.17. 



5 CONCLUDING REMARKS 
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Figure 1: Evolution of the radius for the Schwarzschild type model. 
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Figure 2: Density (multiplied by 10 3 ), energy density flux (multiplied by 10 4 ), 
pressure (multiplied by 10 3 ) and matter velocity (multiplied by 10) for the 
Schwarzschild type model as a function of the time-like coordinate and different 
pieces of the material: 0.2 (solid line); 0.4 (dotted line); 0.6 (small-dashed line); 
0.8 (dashed line); 1.0 (dot-dashed line.) 
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Figure 3: Evolution of the enforced radius A (dashed line) and the resulting 
luminosity E (solid line multiplied by 10 2 ) for the Schwarzschild type model in 
Schwarschild and Bondi coordinates. 



13 



8.00 



< 




5 75 1 1 1 1 1 1 1 

10 20 30 

t 



Figure 4: Evolution of the radius for the Tolman VI type model. 





Figure 5: Density (multiplied by 10 3 ), energy density flux (multiplied by 10 4 ), 
pressure (multiplied by 10 3 ) and matter velocity (multiplied by 10) for the 
Tolman VI type model as a function of the time-like coordinate and different 
pieces of the material: r/a — 0.2 (solid line); 0.4 (dotted line); 0.6 (small-dashed 
line); 0.8 (dashed line); 1.0 (dot-dashed line.) 
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